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ABSTRACT 

Using the relativistic equations of radiation hydrodynamics in the viscous limit, we analyze the 
boundary layers that develop between radiation-dominated jets and their environments. In this paper 
we present the solution for the self-similar, 2-D, plane-parallel two-stream problem, wherein the jet 
and the ambient medium are considered to be separate, interacting fluids, and we compare our results 
to those of previous authors. (In a companion paper we investigate an alternative scenario, known as 
the free-streaming jet model.) Consistent with past findings, we show that the boundary layer that 
develops between the jet and its surroundings creates a region of low-density material. These models 
may be applicable to sources such as super-Eddington tidal disruption events and long gamma-ray 
bursts. 

Subject headings: galaxies: jets — gamma-ray bursts: general — radiation: dynamics — relativistic 
processes 


1. INTRODUCTION 

Astrophysical jets almost certainly exist as aggregates 
of massive particles, magnetic fields, and radiation. In 
certain scenarios, however, the contribution of radiation 
to the energetics of the outflow far outweighs those of 
the particles and magnetic fields, meaning that one can 
essentially neglect the presence of the latter two entities. 

One such scenario occurs in the collapsar model 
of long gamma-ray bursts (GRBs; IWooslevI I1993I 
iMacFadven fc Wooslevl[l999f) . In this model, the core of 
a massive, evolved star collapses directly (or with a short¬ 
lived neutron star phase) to a black hole during the infall 
stage of a type-II supernova. The energy released by the 
material accreting onto the black hole, and ultimately 
observed as the gamma-ray emission, is collimated into 
bipolar jets - the jet formation confirmed by energet¬ 
ics arguments (IWaxman et al.ll 19981 : iFruchter et al.l[r999l : 
iFrail et ^l2001f) and the o bservations of breaks in the X- 

ray a f terglow light curves (lPanaitescull2007t iDado et ^ 
120081 IRacusin et al.l 1200911 - and is often sufficient to 
unbind the stellar envelope, resulting in a supernova 
(IGalama et al.111998I IBersier et al.] 12004 iKamble et al.l 

120091 iLevan et al.112014 but see iFvnbo et al.ll2006ll ! If 

one assumes that the mass of the remnant black hole 
is on the order of a few solar masses, its accretion lu¬ 
minosity exceeds the Eddington limit by roughly ten 
orders of magnitude, meaning that radiation pressure, 
even if the flux is nearly isotropic, is likely an impor¬ 
tant mechanism for driving and sustaining the outflow 
(the fireball model; iRees &: Meszaroi 1199211 . Even if 
the jet is launched by magnetohydrodynamical mech- 
anisms (IBlandford fc Znaield 119771: IBlandford fc Payii^ 
I1982I1 . radiation could still play a prominent role in deter¬ 
mining the dynamics of the jet. Arguments concerning 
the time necessary for the jet to break through the stellar 
eiwelope also seem to disfayor Poynting-doniinated jets 
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(jBromberg et al.ir2015l but see lMundell et al.ll2013ll . 

Jets produced during t idal disruption eyents (TDEs; 
iGiannios &: Metzgei1l2011il - when a star is destroyed by 
the tidal force of a supermassiye black hole - could pro- 
yide another class of radiation-dominated outflow. After 
the star is tidally disrupted, roughly half of the shred¬ 
ded debris remains bound to the black hole and returns 
to the tidal disruption radius. If the black hole has 
a mass less than roughly IO^Mq, that rate of return 
can exceed the Eddington limit of the black hole by or¬ 
ders of magnitude for a signif icant amount of time (on 
the order of days to mon ths; lEyans fc Kochaii^ 119891 
iStrubbe fc Quataertll200^ . Provided that this material 
can rapidly accrete onto the black hole, which is likely 
the case due to the tidal dissipation o f kinetic energy 
(|Kochaneklll994t [Gmllodipn e t_aLl[201^ and rela tivistic 

precession effects (pLee sl ll988l: lEvans fc Kochaneklll989ll . 

the energy released during the accretion process will also 
be supercritical. It was during this supercritical phase 
that the event Swift J 1644+57 was seen to have an asso¬ 
ciate d jetted outflow ( Bu rrows et al.ll20l H iBloom et al.l 
1201 It ICannizzo et al.l 201H IZauderer et a.1.11201lh Ithe 
source Swift J2058+05 may pro vide another examp le of 
a jetted, super-Eddington TDE; ICenko et al.ll2012ll . Al¬ 
though the jet launching mechanism for this event is un¬ 
certain, the magnetic field of the tidally-disrupted star, 
assuming its flux is approximately conserved, is almost 
certainly insufficient to power the outflow. Therefore, 
unless one invokes the exis t ence of a fossil ma gnetic field 
(|Tchekhovskov et al.l 1201-3 : IKellv et al.ll20l4l . the radi¬ 
ation pressure associated with the accretion luminosity 
likely plays some role in powering the jet. At any rate, 
the radiation released during the supercritical accretion 
process affects the dynamics of the collimated outflow 
and contributes substantially to its overall energy and 
momentum. 

Collapsar jets inject a significant amount of energy 
into the overlying stellar envelope as they punch their 
way into the circumstellar medium, creating a pres¬ 
surized “cocoon” of shocked material with which the 
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jet interacts (iM orsonv et al.l l2007t iLazzati et alJ 120071 : 
iLopez-Camara et alJ 120131) ■ and the progenitors them¬ 
selves -jjsunlly_t^en_to_be_JugM^^;evx)lve^ Wolf-Rayet 
stars (jMatznei]i2003l : iWoosley fc HeEfeij|2006ll - are likely 
sustained primarily by radiation pressure. The ini¬ 
tial stages of collapsar jet propagation are therefore 
characterized by the transfer of energy and momen¬ 
tum between two radiation-dominated fluids. Because 
of the high accretion rates and low specific angular 
momentum, the fallback disks generated during the 
super-Eddington phase of TDEs are likely optically 
and geomet r ically thick and radia t ion pressure-supported 
llReej lTO^ iLoeb fc UlmedflOQ^ iCouehlin fc BegelmanI 
l2014all. In the zero-Bernoulli acc retion (ZEBRA) model 
of iCoughlin fc BegelmarJ (|2014af) . for example, the ac¬ 
creting material is inflated into a quasi-spherical envelope 
that surrounds the black hole. The dynamics of the jets 
of supercritical TDEs are therefore also modulated by 
the presence of a radiation pressure-supported, external 
medium. 

Previous authors used the supersonic propagation of 
the outflow to model the interaction of the collapsar 
jet with the overlying en velope as an oblique shock- 
boundary layer st r ucture dBromberer fc Le vin sonI Ip^OOTl : 

I Kohler et al.ll 201 ^ iKohler fc BegelmanI 120121 7 1201411 . In 
these models, the outflow is assumed to consist of a per¬ 
fect, single fluid with a relativistic equation of state, 
and s ome have included the pre sence of magnetic fields 
(e.s.. [Kohler fc Begelmaiil 120121 1. A more realistic pic¬ 
ture, however, is obtained by considering the jet as a 
composite of the massive scatterers present in the out¬ 
flow and the radiation that accompanies it. Indeed, this 
approach constitutes the underlying framework of radi¬ 
ation hydrodynamics. Furthermore, non-ideal, i.e., vis¬ 
cous, effects will tend to “smear” the discontinuities oth¬ 
erwise present in the system, resulting in a more gradual 
transition of the fluid quantities between the jet and the 
environment. 

In the limit that the mean free path of a photon is 
small, radiation acts like an effective viscosity, with a co¬ 
efficient of dynamic viscosity that depends both on the 
radiation pressure and the density of scatterers (see sec¬ 
tion 2 of this paper), and transfers mo mentum and en¬ 
ergy between neigh boring fluid elements (|Weinberg|ll971l : 
ILoeb fc LaoilllQ^ . In the radiation-dominated interac¬ 
tion between collapsar jets and TDE jets and their re¬ 
spective ambient medi a, the effects of rad i ation viscosity 
should be quite large. lArav fc BegelmanI (|1992[ 1 consid¬ 
ered the effects of radiation viscosity on the evolution 
of boundary layers in the two-stream approximation (see 
section 3 of this paper); their treatment, however, was of 
a non-relativistic nature, meaning that their results have 
limited applicability to GRB and TDE jets. 

In this paper, the first of two, we consider the effects of 
radiation viscosity on the propagation of relativistic jets 
in radiation-rich environments. In section 2 we present 
the equations of radiation hydrodynamics in the viscous 
limit. Section 3 applies those equations to the two-stream 
problem, wherein the jet and ambient medium are con¬ 
sidered to be two separate, interacting fluids, and we 
compare our resul t s to t he non-relativistic treatment of 
lArav fc BegelmanI l|1992ll . In section 4 we discuss the re¬ 
sults of the analysis and comment upon the application of 
our models to the jets produced by supercritical TDEs 


such as Swift J1644-I-57, GRBs, and other astrophysi- 
cal sy stems. In a second paper (iGoughlin fc BegelmanI 
l2015h we present an alternate model, the free-streaming 
jet solution, and compare it to the two-stream solution 
presented here. 


2. GOVERNING EQUATIONS 

As mentioned in the introduction, radiative forces be¬ 
have like an effective viscosity in the presence of shear, 
when the change in fluid quantities across the mean free 
path of a photon is small. The precise means by which 
this viscous coupling manifests itself can be determined 
b y analyzing the Boltzm a nn equ ation. 

ICoughlin fc Begelmaij (|2014bll used the general rela¬ 
tivistic Boltzmann equation for Thomson scattering to 
discern the effects of radiation viscosity in the presence 
of both relativistic velocities and gravitational fields (i.e., 
in accelerating reference frames). Instead of reproducing 
their work here, we will simply quote the equations of ra¬ 
diation hydrodynamics for a cold gas (gas pressure much 
less than both the gas rest mass density and the radia¬ 
tion pressure) that result from their analysis (see their 
equation 49): 
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Here the speed of light has been set to one, Greek indices 
range from 0 - 3, p' is the fluid rest-frame density of scat¬ 
terers, e' is the fluid rest-frame radiation energy density, 
K is the scattering opacity (in units of cm^ g~^)) gfj.v is the 
metric of the spacetime, is the covariant derivative, 
is the four-velocity of the flow, and = U^^U'^+gt^'' 

is the projection tensor. The Einstein summation con¬ 
vention has been adopted here, meaning that repeated 
upper and lower indices imply summation. This equa¬ 
tion shows that the coefficient of dynamic viscosity, 77 , 
for an optically-thick, radiation-dominated gas is 


8 e' 


( 2 ) 


which agrees with previ ous findings (iBlandford et al.l 
119851 ILoeb fc Laoilli99^ : note that this specific coef¬ 
ficient is only for the case when the gas and radiation 
interact through Thomson scattering. The gas energy 
equation, which will also be useful for us, can be ob¬ 
tained by contracting equation ([H with I/,^, which we can 
show b ecomes (see equation 50 of iCoughlin fc BegelmanI 
I2ni4hll 
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To close the system, we require that the normalization 
of the four-velocity be upheld and that particle flux be 
conserved: 


= - 1 , 


p'Uf^ 


= 0 . 


(4) 

(5) 


Equations © and © - © constitute six linearly inde¬ 
pendent equations for the six unknowns C7^, e', and p'. 

In addition to the energy density of the radiation, e', 
one can also calculate the number density of photons 
by requiring th at the number flux, be conserved. 
One can show (|Coughlin fc BegelmanI l2014bH that the 
equation V = 0 becomes, in the viscous limit. 


AT'C/M 


= V,. 


^ +N'U°‘y aN' 

P'k, \ 9 6 


( 6 ) 


where N' is the rest-frame number density of photons. 
Once we solve the equations of radiation hydrodynamics 
for the four-velocity of the fluid and the mass density of 
scatterers, we can solve equation © for the number flux 
of photons. 

The goal of the next two sections is to apply equations 
o and © - © to the boundary layers established be¬ 
tween fast-moving jets and their environments. For a 
more thorough discussion of the nature of the equations 
of radiation hydrodynamics in the viscous li mit, we refer 
the reader to ICoughlin fc BegelmanI (l2014bf) . 


3. TWO-STREAM BOUNDARY LAYER 

lArav fc BegelmanI l|1992f) considered the lBlasiu^ (jl908fl 
boundary layer problem, wherein one analyzes the dy¬ 
namics of viscous flow over a semi-infinite, rigid plate, 
with the viscous effects attributed to radiation. They 
showed, however, that the requirement that both velocity 
components vanish on the lower plate, the no-slip con¬ 
dition, results in a divergent boundary layer thickness. 
The authors then examined the more physical scenario 
of two interacting fluids, one moving at some asymptotic 
velocity and the other asymptotically stationary, known 
as the two-stream approximation. In this case the no-slip 
condition no longer applies, and they were able to show 
that the boundary layer thickness remains finite. 


The treatment of lArav &: BegelmanI (1199211 was non- 
relativistic, meaning that their results have limited util¬ 
ity when one considers the boundary layers established 
between GRBs and super-Eddington TDEs and their am¬ 
bient media. The enthalpy of the radiation, which is not 
ignorable in radiation-dominated flows, was also not in¬ 
cluded in their momentum equation. Here we extend 
their analysis to the relativistic regime and we include 
the radiation enthalpy. 

3.1. Basic setup 

We assume that the flow is plane-parallel with no vari¬ 
ation in the a:-direction. The covariant derivatives in 
equations m and © - © can therefore be replaced by 
ordinary partial derivatives. Even though this simpli¬ 
fication significantly reduces their complexity, the most 
compact representation of the equations is still given by 
© and (© - © with d/dx^, so we do not write 

them again here. 

The majority of the motion is along the z-axis, mean¬ 
ing Vz ^ Vy and Vx = 0. At some initial point y = z = 0, 
the “jet,” whose constant, asymptotic (y —>■ oo) veloc¬ 
ity is denoted Vj, encounters the ambient medium, the 
asymptotic (y —>■ —oo) velocity of which is zero. The 
asymptotic densities of the jet and the ambient medium, 
denoted p' and respectively, are both taken to be 
constant. The line y = 0 represents the surface that di¬ 
vides the jet and ambient material, and consistent with 
any boundary layer analysis, we also assume that the 
gradient along y is much greater than that along z, so 
d/dz <c d/dy. 

3.2. Boundary layer equations 

The complexity of equations © and © - © can be 
reduced by introducing the boundary layer thickness Sy 
such that 6 ^ byjbz is a small parameter when 8z is 
chosen to be a typical length scale in z. By keeping 
terms only to lowest order in d, we will recover a set of 
reduced boundary layer equations. 

To determine how b depends on asymptotic fluid 
quantities (it is the reciprocal of the square root of 
the Reynolds number in the classical Blasius boundary 
layer), we compare the lowest-order terms in 6 in the 
gas energy equation to the inviscid terms. By equating 
these terms we are requiring that the viscous heating 
of the radiation contribute a non-negligible increase in 
the entropy of the fluid, but because the inviscid terms 
are proportional to the divergence of the four-velocity 
(see equation [1®, this equality can only be true when 
the gas is compressible, i.e., when the flow velocity is 
supersonic. When the flow becomes very subsonic, the 
energy equation can be ignored as the fluid is essentially 
incompressible. Making the substitutions d/dy ^ ^/by, 
d/dz ^ 1/z, V ^ Vj and p' ^ pg in equation ©, where 
Po = p'j if we are in the jet {y > 0) or p/ if we are in 
the ambient medium (p < 0), we find that the boundary 
layer thickness scales as 


p'^KzVjVj' 

where Tj = (1 — is the Lorentz factor of the jet. 

The boundary layer thickness therefore scales roughly as 
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equations for / and g. It is mathematically convenient, 
however, to dehne a new self-similar variable ^ by 


where r ~ PqK-z is the fluid-frame optical depth 

along 

We can now use our expression for <5 ([71) to keep only 
lowest-order terms in equations ([T]) and (jl]) - ([5]). The 
resulting v = y, i/ = z, gas-energy and continuity bound¬ 
ary layer equations are, respectively. 



( 8 ) 
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e'U>^ 


= 0 , 

(9) 

1 ,917^ 8 e' dUf^dUf^ 

3^ dxi^ 27 p'k dy dy 


C = / 9{a)da, (17) 

Jo 

where d is a dummy variable of integration, and write the 
functions / and g in terms of t his variable. This approac h 
is similar to the one taken by lArav fc Begelmanl (|1992f) . 
With this parametrization, the velocities are 




(18) 


p'TVy = 


2yi V 27 


8 Po Tj Vj 


- f 


(19) 


_9 

dz 



_9 

dy 


p'Tvy 


= 0 . 


( 11 ) 


The first of these demonstrates, as in the non-relativistic 
limit, that the pressure is constant across the boundary 
layer. For the remainder of this section we will assume 
that e'{z) = e'g, i.e., that the radiation energy density is 
independent of z. 


3.3. Self-similar approximation 

The solution to the continuity equation m can be 
obtained by introducing the stream function tjj through 
the definitions 


p'kTvz 


dif 


p'nTvy = 


dtp 


( 12 ) 

(13) 


which manifestly solves the continuity equation; we in¬ 
troduced a factor of k, the opacity, to ensure that ip 
is dimensionless. These relations also demonstrate that 
Vy ^ Svz, which is what we expect: the velocity in the 
direction perpendicular to the majority of the motion is 
reduced by a factor of S. Note, however, that the def¬ 
inition of the stream function now involves the density, 
which is not constant in this analysis. 

As is done in the standard Blasius treatment, we as¬ 
sume that the stream function varies self-similarly as 


\/ (14) 

where 


a = 2//52/=(15) 

is our self-similar variable (we used equation (O for the 
boundary layer thickness) and / is a function to be deter¬ 
mined from equations ([5]) and (flUl) . We will also assume 
that the density varies self-similarly as 

p' = Po9{a), (16) 

where g is a second function. 

One can use these definitions in equations (jH]) and (flUl) 
to derive a set of coupled, nonlinear, ordinary differential 


where a function with a subscript ^ denotes the derivative 
of that function with respect to i.e., /j = df /d£,, /jj = 
d'^f/dfP^, etc. Substituting these relations and the self¬ 
similar scaling for p' , equation (fT6)l . into equations ([9]) 
and (fTUl) . we find the following self-similar equations for 
/ and g: 


1 


2„,2 






„ f 3 , , gHh(? 
90 2 + 


pgf^a^ 

( 20 ) 

( 21 ) 


where we have defined p as the ratio Cg/pg. The 
term —2pff^^/3 and the last term on the left-hand 
side of equation (12011 we re absent in the treatment of 
lArav fc Begelmanl ljl992f) because they did not include 
the enthalpy of the radiation. In the non-relativistic, 
p <C 1 limit, equation (I^U)) reduces to the standard Bla¬ 
sius equation (by rescaling the self-similar variable) and 
q and / decouple from one another, as was found by 
lArav fc Begelmanl (|1992ll . 

Before proceeding further, recall that the two-stream 
problem separates the jet and ambient material into two 
distinct media, the dividing line for our problem chosen 
to be ?/ = 0, and that the asymptotic densities attained 
in these two media are p' and p(j for the jet and ambi¬ 
ent material, respectively. Therefore, in equations (|14l) . 
m and (HU), the parameter pg refers to either p' or 
p(j depending on the sign of a and, hence, The func¬ 
tions / and g are thus piecewise dehned about y = 0, 
as are the self-similar variables a and with solutions 
for ^ > 0 corresponding to jet quantities and those for 
^ < 0 corresponding to the ambient medium. Equations 
(j^Dl) and (!?!]) should therefore be interpreted as a total of 
four equations for four functions: / defined in the jet, / 
defined in the ambient medium, g defined in the jet, and 
g defined in the ambient medium. Note that p can also 
differ between the two media, depending on the asymp¬ 
totic density (but Cg must be continuous across y = 0 
because of equation ([5])). 

Keeping in mind this subtlety of equations (1201) and 
m, we must additionally determine the boundary con¬ 
ditions satisfied by / and g. The first two conditions sat¬ 
isfied by / can be determined by recalling that v —>■ Vj as 
y —> 00 and n— >0 as ?/—>■—oo. Investigating equation 
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m and noting that ^ scales with y, these requirements 
translate to 


/ 4 (oo) = l, / 5 (-oo) = 0. (22) 

We also require that the density approach its asymptotic 
values in the limits of y —?> ±oo. From equation (fTHl) . this 
gives 


g(±oo) = 1. (23) 

Now, note that if the jet and ambient materials are to 
retain their respective identities, then the flow along the 
surface of contact at j/ = 0 must remain parallel to that 
surface. In other words, there must not be any mass flow 
across the boundary, i.e., this surface is a contact dis¬ 
continuity, which means that ny(0) = 0. From equation 
m, this shows that / must satisfy 


/(O) = 0. (24) 

The other two boundary conditions can be determined 
by requiring that the normal and transverse components 
of the energy-momentum tensor, the divergence of which 
we took to obtain equation o, be continuous across the 
point of contact y = 0. We can show that these restric¬ 
tions demand that and be continuous across y = 0, 
which closes the system. 

We can also calculate the comoving-frame number den¬ 
sity of photons N' by solving equation (IH]). Keeping only 
terms to lowest order in the boundary layer thickness, 
the equation of conservation of photon number becomes 




N'Uf^ 


1 d / I dN'\ 
3 dy \p'k dy ) 


= 0 . 


(25) 


As we did for the number density of scatterers, we will as¬ 
sume that the photon number density varies self-similarly 
across the boundary layer as 


N' = N'h{^), (26) 

where Nq is the asymptotic number density of the jet 
or the ambient medium. With this form for the number 
density, we find that the equation for h becomes 


3p2 2 

2 n + 


h-fh 



(27) 


As for the functions / and g, h is really piecewise defined 
across the boundary j/ = 0, and so equation (1271) is really 
two equations - one for the photon number density in the 
jet, and another for that in the ambient medium. The 
requirement that the number density of photons asymp¬ 
totically approach its jet and ambient values gives 


to the physical coordinates y and z (equation (fT51) l. The 
transformation can be achieved by differentiating equa¬ 
tion m, rearranging, and integrating to yield 



= a, 


(29) 


where ^ is an integration variable and we set the integra¬ 
tion constant to zero because we demand ^(a = 0) = 0. 
After solving equations (pOl) and ((211) for g{^), we can nu¬ 
merically integrate and solve equation (1^^ for ^ in terms 
of a. 


3.4. Solutions 

Equations (Uni), dH]) and (EH), together with the 
boundary conditions on /, g, and h, govern the behavior 
of the velocity, number density of scatterers, and number 
density of photons throughout the two-stream boundary 
layer. Here we present and analyze the solutions to those 
equations as we vary the asymptotic jet Lorentz factor 
Fj and the quantity p = Bq/pq. 

One caveat with the definition of a is that it depends 
on Fj via equation ((Tbl) . Therefore, if we plot solutions 
with different F^, we must be careful to incorporate this 
dependence so that the range of physical space we con¬ 
sider for each solution is the same. Because of this fact, 
in the following figures we will plot our solutions as func¬ 
tions of the variable 


a 

a = 


\/i V 8 


V . 27 


PoK- 


(30) 



a 

Fig. 1.— The function /^, which is the normalized ^-component 
of the four velocity, in terms of the parameter d oc yl^/z, for /i = 1 
and a number of jet Lorentz factors, as indicated by the legend. As 
we can see, the thickness of the velocity boundary layer, in terms 
of d, is nearly independent of Tj. 


h{±oo) = 1. (28) 

The y-component of the flux must also be continuous 
across the boundary, which, when written out, shows 
that the derivative of h must be continuous across ^ = 0, 
which yields the final two boundary conditions. 

Finally, the solutions to equations (1^ . (1^ . and (l?7)) 
will be in terms of the parameter which is itself a func¬ 
tion of g (equation (fTTll l. We would like the solutions to 
be in terms of the parameter a, which is directly related 


Figure [T] shows the solution for the normalized z- 
component of the four-velocity (/^) as a function of the 
self-similar variable a oc yjy/z for a number of jet Lorentz 
factors and p = 1. We see that the width of the boundary 
layer is nearly unchanged as we vary the Lorentz factor 
of the asymptotic jet. The value of /j at the contact 
discontinuity is lower for larger F^-, resulting in a greater 
shear (^ /jj) as one proceeds into the jet. The flatten¬ 
ing of the velocity around 5 = 0 arises from the behavior 
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tion shows that g oc — l/lii(^). 

In Figure 13] we plot the solution to equation (l?7l) . the 
normalized, rest-frame number density of photons, for 
the same set of Lorentz factors and p = 1. By compar¬ 
ing this with Figure 121 we see that the photon number 
density roughly tracks that of the scatterers. However, 
the photon number density remains below its asymptotic 
value significantly farther into the ambient medium than 
does the particle density. Because the energy density of 
the radiation is constant across the boundary layer, the 
decrease in photon density corresponds to a higher av¬ 
erage energy per photon increases within the boundary 
la yer - a clear manifestation of viscous heating, as noted 
by lArav &: BegelmanI l|1992ll . 


Fig. 2.— The variation of the normalized density, given by g, 
as a function of a for the same parameters as those chosen for 
Figure^ The density remains below its asymptotic value over a 
slightly larger range of a for higher Fj, and the decrease in den- 
sity within the boundary layer is consistent with the findings of 
lArav &: BegelmanI <119921) . The density formally equals zero at 
d = 0; however, because g approaches zero at a very slow rate 
(recall g{^) oc —1/ln^), it appears from the Figure, which only 
samples a finite number of points around d = 0, that the density 
remains positive and larger for smaller . 



a 

Fig. 3.— The solution to equation dSTJi, h, which is the normal¬ 
ized number density of photons, for the same set of parameters as 
those chosen in Figure ^ The number density of photons is seen 
to roughly track the number density of scatterers. Because the 
energy density of the radiation remains constant across the layer, 
the energy per photon increases in the boundary layer. 


of the density around this region and the function ^(a) 
determined therefrom. 

Figure [5] illustrates the manner in which the density 
varies over the boundary layer for the same parameters 
as those chosen for Fi g ure [T1 Consistent with the findings 
of lArav fc BegelmanI (jl99^ . we find that the transition 
from the jet to the ambient medium carves out a region 
of low density material. This behavior can be under¬ 
stood by noting that the shear in the flow causes viscous 
heating of the fluid, which results in an increase in the 
specific entropy s'. Since the specific entropy scales as 
s' oc for a radiation-dominated gas, an increase 

in the entropy at constant pressure corresponds to an 
decrease in the density of scatterers. As was found by 
lArav fc Begelma^ (|1992f ). the density equals zero at the 
origin, which can be gleaned from the asymptotic behav¬ 
ior of equation (1^ : for ^ <C f, we can let ^ ^ 5, 

and solving the resultant approximate differential equa¬ 



a 

Fig. 4. — The solution for the normalized ^-component of the 
four-velocity (/^) for Fj = 25 and three different values of /i., as 
indicated by the legend. The velocity profile does not differ much 
as n increases beyond 1, but for smaller values of fi the boundary 
layer widens noticeably. 



a 

Fig. 5.— The function g, which is the rest-frame number density 
of scatterers, for the same parameters as in Figure [4| For smaller 
values of /i, the density is significantly reduced from its asymptotic 
value over a larger range in a. 


Figures 131151 and [51 show the z-component of the four- 
velocity, the number density of scatterers, and the photon 
number density, respectively, for Tj = 25 and /i = O.f, 1, 
and 10. The overall qualitative behavior of the solutions 
is similar to that depicted in Figures [Tl-|31 More specif¬ 
ically, however, we find that values of p greater than 
unity cause the boundary layer thickness to decrease, 
but not appreciably. On the contrary, a value of p only 
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a 

Fig. 6.— The solution to equation m, which gives the number 
density of photons observed in the comoving frame, for the same 
set of parameters as in Figure |4] The radiation number density 
roughly follows that of the scatterers. 

marginally less than one results in a significant widening 
of the boundary layer. This dependence is ultimately re¬ 
lated to the compressibility of the fluid and the relation 
between that compressibility and the sound speed of a 
radiation-dominated gas (see discussion below). 

4. DISCUSSION 

The plots of the previous subsection demonstrate how 
radiation viscosity affects the boundary between a fast- 
moving flow and its surroundings under the two-stream 
approximation. Our anal ysis generalizes the treatment of 
lArav fc BegelmanI (|1992l l by permitting the jet velocity 
to be relativistic and by incorporating the enthalpy of 
the radiation in the moment um equation. Our results 
are similar to those found by lArav fc BegelmanI (Il992f) 
(compare, e.g., their Figure 3 to our Figures [T] and [J); 
there are, however, a few notable differences. 

For one, our boundary laye r thickness, given by equa - 
tion ©, differs from that of lArav fc Beedmanl l|1992f) . 
who found S'^ ^ e'/{p'o^KVjz). In addition to the Lorentz 
factor contained in ours, their boundary layer thickness 
has an additional factor of e'jpQ, meaning that, in the 
non-relativistic limit, our results do not agree. This dis¬ 
crepancy arises from the fact that, while we compared the 
lowest-order terms in 6 to the inviscid terms in the gas 
energy equation to obtain our boundary layer thickness, 
they compared the viscous term to the inertial term - the 
one proportional to p' - in the z-component of the mo¬ 
mentum equation to obtain theirs. Because they ignored 
the enthalpy of the radiation, the inertial term was the 
only inviscid term present in the momentum equation, 
making it the only one available to balance the viscous 
part. However, if one does not ignore the radiative contri¬ 
bution to the momentum equation, one can now equate 
the viscous terms to either the term proportional to the 
mass density or the one proportional to the enthalpy. 
Because the density equals zero at a = 0, there is al¬ 
ways some location at which the enthalpy term exceeds 
the inertial term, making it more reasonable to equate 
the former to the viscous contribution than the latter. If 
one follows this route, one recovers our ordering for the 
boundary layer thickness. 

The second difference is that their solutions depend 
on the square of the Mach number, = v'^jcg, where 
Cg = 4e'/ (9/9q) is the non-relativistic sound speed. Our 


solutions, on the other hand, depend on both the jet 
velocity Vj and the ratio p = e'jp^, which is propor¬ 
tional to the non-relativistic sound speed. One reason 
for this difference arises from the discrepancy between 
our definitions of the boundary layer thickness. Another 
is due to the fact that our solutions included the enthalpy 
of the radiation; had lArav &: Begelma^ (Il992f) included 
this term, factors of the sound speed would have arisen in 
their z-momentum equation. Finally, our solutions also 
depend on the value of the jet velocity - not just the 
ratio of the jet speed to the sound speed - because we 
included all relativistic effects, meaning that the speed 
of light now plays a role in determining the evolution of 
the system. 

Our solutions show that, for fixed p, the thickness of 
the boundary layer is approximately independent of the 
jet Lorentz factor, which is due to the competition be¬ 
tween the scaling of the fiducial boundary layer thick¬ 
ness 6 oc \/ (see equation (O) and the viscous heat¬ 
ing. Specifically, a larger Fj results in a smaller 5 and 
a greater shear; this shear increases the specific entropy 
s' oc e'/(p')'*/^ and, since the pressure is a constant, this 
increase in the entropy implies a decrease in the density 
of scatterers which widens the boundary layer. 

Figures m [5l and [6] show how the solutions for the ve¬ 
locity, density of scatterers, and density of photons vary 
for a fixed jet Lorentz factor (we chose F^ = 25) but for a 
variable p. Increasing p relative to p = 1 tends to slightly 
decrease the boundary layer thickness, while decreasing 
the value of p relative to p = 1 dramatically increases 
the thickness. This behavior arises from the fact that 
the viscous heating, which decreases the density of scat¬ 
terers and widens the boundary layer, is effective when 
the gas is compressible. The compressibility of the fluid, 
however, is only important when the flow velocity is su- 
person ic, and we can show that the sound speed of the 
gas is (jCoughlin fc Begelmanll^l4b|l 


-2 / M 
3Vl + 4p/3’ 


(31) 


which, noting that e' — 3p', where p' is the radiation 
pressure, correctly reduces to Cg — V^4p'/(3p') in the 
limit p 1 and Cs ~ I/a/S in the limit p ^ 1. Thus, 
when p ^ 1, the location of the sonic point extends 
farther into the ambient medium, resulting in a widening 
of the boundary layer. On the other hand, when p ^ 1, 
the sonic point approaches the jet, but only slightly due 
to the fact that the sound speed approaches a constant 
as p —>■ c». In fact, based on this reasoning, we expect 
our solutions to be independent of p in the large p limit, 
which is indeed reflected in equation (l20ll . 

The densities of scatterers and photons both decrease 
within the boundary layer. This behavior has two inter¬ 
esting consequences. First, the lower number density of 
scatterers means that the optical depth is lower in the 
boundary layer. We are therefore able to see farther into 
the medium along lines of sight that probe regions of 
high shear. Second, because the energy density of the 
radiation stays unchanged as we traverse the media, the 
average energy per photon increases, resulting in a harder 
spectrum within the boundary layer. 

Equation ([8]) shows that the radiation energy density, 

























and hence the pressure, is constant across the bound¬ 
ary layer, which is ultimately a statement of the causal 
connectedness of the jet. This means, equivalently, that 
the boundary layer thickness 5y can be traversed by a 
sound wave in less time than it takes the jet material to 
cross the distance 5z. Since the transverse sound speed 
is related to the isotropic sound speed by Cs_l = c-s/Tj, 
we find that the boundary layer thickness 5 must satisfy 
5 < Cs/Tj. When this inequality is no longer satisfied, 
equation ([8|) no longer holds, and we must include more 
terms in all of the boundary layer equations that account 
for changes in the pressure. 

The solution for g equals zero at a = 0 which, as we 
mentioned, can be determined by inspecting the asymp¬ 
totic —>■ 0) limit of equation (1211). This feature was 
also found by lArav fc Besrelmanl (Il992| l , and can be un¬ 
derstood physically by noting that, when the jet and the 
ambient medium initially interacted, the surface separat¬ 
ing them served as a discontinuity in velocity, resulting 
in an infinite shear and entropy generation. Although 
the solutions presented here illustrate the time-steady 
state of the system after radiation viscosity has smoothed 
out the discontinuity, the infinite entropy along the con¬ 
tact discontinuity persists and drives the density to zero. 
Mathematically, this behavior is ultimately due to the 
fact that /(O) = 0, which itself came from the require¬ 
ment that there be no mixing at the surface separating 
the jet and the external environment. This boundary 
condition is necessary to ensure that the two media retain 
their respective identities and underlies the two-stream 
assumption, and it allows us to prescribe different asymp¬ 
totic properties of those media. 

5. SUMMARY AND CONCLUSIONS 

We applied the relativistic equations of radiation hy¬ 
drodynamics in the viscous limit to the two-stream 
boundary layer, expan ding on the past work of 
lArav fc Begelmaa (Il992f) . These equations are applica¬ 
ble as long as changes in the fluid quantities are small 
over the mean free path of a photon. 

An interesting feature of the solutions presented here 
is the depression in the number density of scatterers 
within the boundary layer separating the jet and the 
a mbient medium, wh i ch is consistent with the findings 
of lArav fc BegelmanI l|1992f) . We also showed that the 
number density of photons N' roughly tracks the density 
of scatterers, reaching a minimum towards, but not at, 
the contact discontinuity. Therefore, observers viewing 
a relativistic outflow with lines of sight that probe re¬ 
gions of high shear see farther into the outflow and they 
also see a higher energy per photon, and hence a harder 
spectrum. 

Our solutions show that the physical boundary layer 
thickness does not depend strongly on the jet Lorentz 
factor Tj, which results from a competition between vis¬ 
cous heating and the scaling of the fiducial boundary 
layer thickness (5^1/ The dependence of our solu¬ 
tions on fi arises from the fact that the change in entropy 
of the flow is related to its compressibility, which is most 
important where the outflow velocity is supersonic. Since 
the sound speed scales as Cg y/Ji when p ^ 1, the point 
at which the outflow becomes subsonic extends farther 
into the ambient medium when /i is small, resulting in a 
widening of the boundary layer. When p becomes larger 


than one, however, the sound speed does not increase 
much, asymptotically approaching 1 /v^, meaning that 
the sonic point only slightly approaches the jet, yielding 
a marginal decrease in the boundary layer thickness. 

A number of assumptions about the nature of the jet 
and its surroundings are built in to the two-stream so¬ 
lutions presented here. Specifically, we adopted a two- 
dimensional, plane-parallel geometry for the flow and its 
surroundings, and we demanded that there be no pres¬ 
sure gradient (e'(z) = Cq) in the ambient medium. While 
these assumptions greatly enhanced the tractability of 
the problem, they somewhat hinder the astrophysical rel¬ 
evance of the solutions, as no systems likely conform pre¬ 
cisely to these restrictions. However, in local regions of 
an outflow, where the radius of curvature is large and the 
pressure can be considered relatively constant, these so¬ 
lutions may be actualized. Therefore, while the physics 
of an entire jet-disk system may be poorly modeled by the 
two-stream scenario, local patches of the outflow, where 
the geometry can be treated as locally flat and the pres¬ 
sure gradient is small, are likely well-described by the 
solutions presented here. 

The two-stream solutions could be applied to rela¬ 
tivistic, radiation dominated jets, such as those that 
appear during super-Eddington TDEs, the event Swift 
J1644-I-57 being one such case. The event Swift 
J2058+05, observed shortly after Swift J1644, is an- 
other source that h as been interpreted as a jetted TDE 
(|Cenko et al.ll20l3l . A comparison between the models 
presented here and these sources could provide valuable 
information concerning their progenitors and the prop¬ 
erties of the jets themselves. The application of these 
models to long GRBs could likewise prove fruitful, po¬ 
tentially yielding, for example, information concerning 
the Lorentz factor of the jet, the stellar progenitor, and 
the direction of the line of sight between the observer and 
the source. 

These models may also be appl icable to jetted X-ray 
binar y systems, or “m i croqu as ars” (Mirabel fc Rodrfgued 
H)99^_jFenderjet alj__ [2004 12009 1. For example. 


Arav fc BeefelmanT i 199,*^ a pplied their non-relativi stic, 


radiation-viscous solution (Arav fc: Begelni^ Il99 2f) t o 
the source SS 433 (jFabrikall2004 iBegelman et al.ll2006^ . 
Since the jets of SS 43 3 are mildly relativistic, their 
speeds being Vj ~ 0.26 c (iMargon fc Andersonlll989ll . in¬ 
cluding the relativistic terms arising from the treatment 
presented here may place new constraints on the proper¬ 
ties of those jets and the surrounding envelope. 

Finally, quasi-stars - protogalactic gas clouds sup- 
port ed by an accreting black hole - may also contai n 
jets (|Begelman et al.l feoOGbl [^008t iCzernv et al.l I2012H . 
Since the power radiated by the black hole at the center 
of a quasi-star supports the overlying gaseous envelope, 
the mass of the envelope greatly exceeding that of the 
hole, the black hole accretes supercritically by several or¬ 
ders of magnitude. The jets launched from these systems 
are therefore radiation-dominated, and as they propa¬ 
gate through the quasi-star envelope, radiation-viscous 
effects likely dominate the interaction between the two 
media. Applying the solutions presented here to these 
systems could then yield information about the proper¬ 
ties of these jets and the role they may have played during 
the epoch of galaxy formation. 
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One drawback to these models and their physical ap¬ 
plication, however, is that they require that the jet and 
ambient medium be separated by a contact discontinuity, 
which results in the non-physical vanishing of the density 
of scatterers at the interface. Furthermore, this prevents 
the jet from entraining ambient material; while this is not 
particularly problematic for the two-stream problem, in 
which the jet is considered infinite in extent, realistic jets 
have a finite width and total momentum, meaning that 
the entrainment of ambient material will cause a decrease 
in the outflow velocity with z that cannot be captured 

with the two-stream treat ment. _ 

In a companion paper (|Coughlin fc BegelmanI [201511 . 


we investigate a different type of viscous boundary layer 
- the free-streaming jet model. This model treats the 
entire system, jet and ambient medium, as a single fluid, 
which removes the need for a contact discontinuity and 
allows the density to remain non-zero throughout the 
boundary layer. We also show that the entrainment of 
ambient material causes an overall slowing of the outflow. 


This work was supported in part by NASA Astro¬ 
physics Theory Program grant NNX14AB37G, NSF 
grant AST-1411879, and NASA’s Fermi Guest Investi¬ 
gator Program. 
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